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PREFACE 


This  Report  is  part  of  Rand's  study  for  the  Advanced  Research 
Projects  Agency  of  those  phenomena  which  affect  the  performance  of 
underwater  optical  reconnaissance  and  guidance  equipment.  The  ob- 
jective  of  these  studies  is  to  provide  sufficient  understanding  to 
permit  the  systems  analyst  to  compute  performance  estimates  under 
various  operational  conditions. 

A  better  understanding  of  underwater  visibility  will  be  obtained 
only  through  an  investigation  of  those  mechanisms  which  give  rise  to 
the  very  intense  small-angle  forward  scattering  of  light  by  ocean 
water.  A  recent  Rand  study  (Ref.  1)  of  the  scattering  by  suspended 
particles  was  incorrect  in  that  it  confused  Individual  (incoherent) 
and  collective  (coherent)  scattering.  These  mechanisms  differ  con¬ 
siderably  in  their  dependence  on  the  particle  density,  and  the  use 
of  the  expressions  given  in  Ref.  1  would  tend  to  underestimate  the 
beam  spreading  for  typical  particle  densities  by  several  orders  of 
magnitude. 

The  present  Report  attempts  to  clarify  some  of  the  distinctions 
between  various  scattering  mechanisms,  and  should  be  of  interest  to 
those  concerned  with  underwater  visibility  and  the  use  of  lasers  in 
underwater  detection. 

This  material  has  been  submitted  to  the  Applied  Optics  Journal 
as  a  correction  to  Ref.  1,  which  appeared  in  the  January  1.971  issue. 

A  revised  version  of  Ref.  1  is  currently  in  process. 


COHERENT  AND  INCOHERENT  MULTIPLE  SCATTERING  OF  LIGHT  IN  SEA  WATER 


Scattering  of  light  in  sea  water  is  dominated  by  scattering  from 
suspended  particles  and  biological  organisms.  Since  the  wavelength  of 
the  radiation  is  small  compared  to  the  mean  distance  between  scat¬ 
tering  centers,  the  particles  scatter  independently  (incoherent 
multiple  scattering).  Coherent  scattering,  important  only  when  the 
wavelength  is  large  compared  to  the  mean  separation,  is  negligible 
for  scattering  by  sea  suspensions.  A  recent  discussion' **'  of  the  theory 
of  multiple  scattering  of  light  by  suspended  particles  is  incorrect  in 
that  it  uses  a  result  from  coherent  scattering  theory  to  calculate  the 
incoherent  scattering  and,  as  a  consequence,  underestimates  the  scat¬ 
tering  in  sea  water  by  several  orders  of  magnitude . 

Both  coherent  and  Incoherent  multiple  scattering  have  been  treated 

it 

exhaustively  in  the  literature  as  they  apply  to  neutron  diffusion, 
cosmic-ray  showers,  stellar  radiative  equilibrium,  ohmic  resistivity, 
electron  penetration,  pulaar  scintillation,  etc.,  and  most  of  the 
general  underlying  theory  is  readily  adaptable  to  scattering  of  light 
by  sea  water.  However,  the  physical  conditions  for  sea-water  scattering 
are  highly  variable,  and  there  is  little  value  in  extensive  treatments; 
at  present,  only  the  simpler  concepts  can  be  experimentally  tested. 
Although  these  concepts  are  generally  well  known  and  have  been  long 
applied  to  the  subjects  mentioned  above,  the  recent  publication  of  an 
Incorrect  treatment^  suggests  that  it  may  be  useful  at  this  time  to 
review  the  simpler  aspects  of  multiple  scattering  theory. 

Notation  will  differ  somewhat  from  that  of  Ref.  1,  but  the  model 
is  the  same:  a  dilute  dispersion  of  spherical  (or  nearly  spherical) 
transparent  scatterers  of  radius  a  and  index  of  refraction  n.  It  is 
assumed  that  ka  »  1,  where  k  is  the  wave  number  of  the  incident 
light.  N(r)  represents  the  number  of  particles  per  unit  volume  at 


Incoherent  single -particle  scattering  is  sometimes  called 
microscopic  or  individual  scattering;  coherent  collective  scattering 
is  sometimes  called  macroscopic  or  bulk  scattering. 
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r;  the  fractional  volume  f(r)  occupied  by  the  spheres  is 

f(0  -  j-  a3  Hr)  (1) 

Incoherent  scattering  refers  to  scattering  by  individual 
particles;  coherent  scattering  refers  to  scattering  by  clusters 
of  particles.  In  either  case,  if  (0*)  represents  the  mean  square 
scattering  angle  from  a  single  particle  (or  a  single  cluster),  then, 
if  the  scatterings  are  independent  events,  the  mean  square  scattering 
angle  (6  )  at  z  is 


<02>  -  p<02>  (2) 

where  p(z)  is  the  mean  number  of  scatterings  associated  with  the 
parameter  z.  If  z  represents  the  penetration  distance,  p  may  be 
expressed  in  the  form 


p  -  <N>  or  (3) 

where  (N)  is  the  mean  particle  (or  cluster)  density  and  o  is  the 

single-particle  (or  ■ ingle- cluster)  total  scattering  cross  section. 

2 

Thus,  in  general,  (0  )  at  z  is  given  by 

<02>  -  <N>  cz  <02>  (4) 

(Higher-order  coherent  effects  may  be  generated  by  clusters  of  clusters, 
etc.)  Henceforth,  in  order  to  distinguish  between  particle  and  cluster 
scattering,  th*.  subscript  c  will  be  used  to  designate  cluster  parameters. 


Equation  (2)  is  the  small-angle  approximation  to  the  more  exact 
.-elation  given  by  (cos  0)  ■  (»!  cos  0^),  where  the  product  extends  over 
all  collisions  and  0^  is  the  i'th  single  scattering. 
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The  parcicle8  are  aisumed  to  be  randomly  distributed  with  density 

N(r)  -  <N>  +  6N(r)  (5) 

where  the  mean  (N)  is  constant  and 

<6N(1)  6N(2)>  -  <(6N)2>  B(r)  (6) 

where  <(6N)2>  is  constant,  r  ■  -  r2l,  and  B(0)  -  1  (the  notation 

fiN(l)  stands  for  6N(r^),  etc.)*  The  characteristic  cluster  size  is 
given  by  the  correlation  length  associated  with  the  two-point  cor¬ 
relation  function  B(r).  Similarly,  f(r)  may  be  expressed  in  the  form 

f(r)  -  <f>  +  6f(r)  (7) 

where,  from  Eqs.  (1)  and  (6), 

<6f(l)  6f<2)>  -  ««f)2>  B(r)  (8) 

The  coherent  properties  of  a  suspension  may  be  characterized  by 
a  coherent  refractive  index  n£(r)  given  by 

n  (r)  -  n  +  u  (r)  ■  n  +  <u  >  +  (r)  (9) 

c  o  c  o  c  c 


where  nQ  is  the  index  (assumed  constant)  when  N  ■  0,  and  where  ^c(r) 
is  given,  for  sufficiently  dilute  solutions,  by 


ft 

The  relation  expressed  in  Eq.  (10)  is  a  familiar  result  applicable 
to  small  (ka  «  1)  dipole  scatterers;  its  use  where  the  particle  scat¬ 
tering  is  predominantly  forward  and  the  particles  are  large  (ka  »  1) 
requires  some  justification.  It  can  be  shown,  for  example,  that  when¬ 
ever  the  (particle)  backscatterlng  is  small  compared  to  the  forward 
scattering,  the  suspension  exhibits  a  coherent  magnetic  permeability 
equal  to  its  coherent  electric  permittivity  ng  (even  though  the  indi¬ 
vidual  particles  are  nonmagnetic).  Only  if  |^ka|  «  1  is  Eq.  (10) 
valid;  under  these  conditions  the  coherent  scattering  cross  section  is 
four  times  that  determined  by  fluctuations  in  permittivity  alone.  For 
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uc(r)  -  n,f(r)  (10) 

where  m>»  given  by 

H  -  n  -  nQ  (11) 

is  a  constant  and  represents  the  difference  between  the  refractive 
index  of  the  scattering  particle  and  that  of  the  medium  (here  the  sea) 
in  the  absence  of  suspended  particles. 

The  mean  coherent  index  of  refraction  of  the  suspension  Is  nQ  +  (p,c) 
the  coherent  scattering  (for  angles  6  ^  0)  is  determined  by  the  fluctua- 
ting  component  6^(r).  FromEqs.  (8)  and  (10) 

<6nc(l)  &UC<2)>  -  <(&nc)2>  B(r)  (12) 

where,  from  Eq.  (1) , 

«S)2>  “  (^*3>A)2  ((6N)2>  (13) 

(Of  course,  other  effects,  temperature  and  salinity  variations,  for 
example,  can  cause  large-scale  fluctuations,  or  patches,  in  the  index 
(i.c.,  nQ  may  not  be  constant).  Light  scattering  by  patches  is 
determined  in  terns  of  the  patch  correlation  function  in  a  manner 
Identical  to  that  used  to  determine  cluster  scattering  from  the  cluster 
correlation  function.) 

Henceforth,  in  the  interest  of  conciseness,  numerical  constants 
of  order  unity  will  be  discarded;  the  mathematics  needed  to  generate 


most  sea  particles,  |uka|  -»  1,  and  Eq.  (10)  is  invalid  (for  example, 

,c  is  then  imaginary  and  independent  of  u);  nevertheless,  the  coherent 
scattering  remains  negligible  compared  to  the  incoherent  scattering 
whenever  the  inter-particle  distance  Is  large  compared  to  the  wave¬ 
length.  To  avoid  a  lengthy  discussion  of  various  aspects  of  the  always 
negligible  coherent  scattering,  and  because  the  relation  expressed  in 
Eq.  (10)  is  adopted  without  comment  in  Ref.  1  (see  Eq.  (7)  of  Ref.  1), 
the  coherent  scattering  is  evaluated  from  Eq.  (10). 


the  missing  constants,  although  straightforward,  is  often  lengthy  and 

is  generally  available  in  the  scattering  literature.  Thus,  in  order 

to  estimate  (8  )  from  Eq.  (4)  using  the  Born  approximation,  the  most 

convenient  approach  is  to  use  the  equivalence  (except  for  numerical 

factors  of  order  unity)  of  the  well-known  geometrical  approximation 

and  the  Born  approximation  in  the  evaluation  of  the  product  of  c  and 
2 

(8.).  For  scattering  by  spheres  of  radius  a  and  index  deviation  ia, 

•  2  2  2 
the  geometrical  approximation  a  ■  Tta  ,  (6^)  ~  p,  ,  gives,  from  Eq.  (4), 

<82)  -  <N>  a2zn2  ~  (f>  (i2  J  (14) 

2  “2 

(It  may  be  shown  that,  in  the  Born  approximation,  (8,)  ~  (ka)  and 

2  2  4  A 

a  ~  pi  k  a  ,  and  Eq.  (14)  is  again  obtained.) 

Equation  (14)  gives  the  Incoherent  scattering.  The  coherent 

scattering  is  determined  by  the  cluster  size  R  (the  correlation  length) 

.3 

and  the  index  deviation  6p  .  The  cluster  density  is  N  ~  R  ,  and,  in 

c  2  2  c  2 

the  geometrical  approximation,  o  ~  R  and  (8,)  ~  ((6p  )  ).  Thus, 

C  a  C  C 

from  Eq.  (4) , 

<e2)c  -  <(6nc)2)  J  (15) 

(In  the  Born  approximation,  (82)£  ~  (l®)  2  end  ac  ~  ((6pc)2)  k2R^.) 
According  to  Eq.  (10),  (82)c  ®*y  be  expressed  In  the  form 

<82>c  ~((6f)2)  p2  2.  (16) 

If  p  represents  the  ratio  of  the  coherent  to  Incoherent  scattering, 
that  is, 
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J 
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Mi 

> 

■ "  T  -*v 


then,  fromEqs.  (14)  and  (16), 
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p  ~  a(f>  £ 


where  or  represents  the  squared  relative  particle  density  fluctuation: 

a  “ 

Since  R  >  d,  where  d  represents  the  average  separation  between  the 
★ 

particles, 


<<W2> 

<N>2 


"  <«*  d~“(d)  <17> 

-8  „ii 

Typically,  f  10  ,  and  hence  p  £  10  for  100  percent  fluctuation; 

that  Is,  the  coherent  scattering  Is  completely  negligible. 

The  previous  exposition^  attempted  to  determine  the  Incoherent 

scattering  by  using  the  formulas  for  coherent  scattering  with  a  cor- 

2 

relation  function  whose  strength  was  Pc ,  with  the  correlation  distance 
given  by  the  radius  of  the  particle.  This  procedure  Is  physically 
meaningless  and  gives  neither  the  incoherent  nor  the  coherent  scat¬ 
tering.  It  was  assumed ^  that  p.  was  the  random  variable  and  that  f 
was  constant;  hence  the  result  is  the  coherent  expression  given  by 
Eq.  (15)  with  6p.c  replaced  by  ^  (■  ^f)  and  R  replaced  by  a: 

(0  )  —  p.  f  z/a 

This  result  (Eq.  (14)  of  Ref.  1)  is  considerably  smaller  (by  the 


Equation  (17)  is  the  large -particle  (ka  »  1)  approximation  to 
the  more  general  result,  valid  for  particles  of  arbitrary  size,  given 
by 


P 


.  2 

a 

~  O'.  — 


♦  x-2'. 
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factor  f)  than  that  (Eq.  (14)  above)  giving  the  dominant  scattering 
* 

process . 


Substitution  of  the  parameters  and  distances  given  in  Ref.  1  into 
the  incoherent  multiple  scattering  formula  given  above  (Eq.  (14)) 
Implies  that  (0^)  »  1.  As  pointed  out  in  the  footnote  on  page  2, 

Eq.  (14)  is  a  small-angle  approximation,  and  hence  ia  applicable  only 
to  distances  for  which  (9^)  ~  1;  at  greater  distances  a  diffusion 
analysis  is  more  appropriate.  Also,  since  Ref.  1  uses  the  incorrect 
scattering  expression  to  determine  f  from  experimental  values  of 
(0Z),  the  values  of  f  given  there  are  too  large  by  several  orders  of 
magnitude.  [R.  F.  Lutomlrakl  (private  communication)  has  pointed  out 
that  the  correct  dependence  on  f  for  the  incoherent  scattering  can  be 
obtained  from  expressions  for  coherent  scattering  with  correlation 
functions  of  the  form  (6p,(l)6u(2))  ■  ^fB(r).] 


« 
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